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We construct static soliton solutions with non-zero Hopf topological charges to a theory which is 
an extension of the Skyrme-Faddeev model by the addition of a further quartic term in derivatives. 
We use an axially symmetric ansatz based on toroidal coordinates, and solve the resulting two 
coupled non-linear partial differential equations in two variables by a successive over-relaxation 
(SOR) method. We construct numerical solutions with Hopf charge up to four, and calculate their 
analytical behavior in some limiting cases. The solutions present an interesting behavior under the 
changes of a special combination of the coupling constants of the quartic terms. Their energies 
and sizes tend to zero as that combination approaches a particular special value. We calculate the 
equivalent of the Vakulenko and Kapitanskii energy bound for the theory and find that it vanishes 
^ I at that same special value of the coupling constants. In addition, the model presents an integrable 

. sector with an infinite number of local conserved currents which apparently are not related to 

symmetries of the action. In the intersection of those two special sectors the theory possesses exact 
vortex solutions (static and time dependent) which were constructed in a previous paper by one of 
the authors. It is believed that such model describes some aspects of the low energy limit of the 
pure SU(2) Yang- Mills theory, and our results may be important in identifying important structures 
in that strong coupling regime. 
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I. INTRODUCTION 



. We construct static soliton solutions, carrying non-trivial Hopf topological charges, for a field theory that has found 
^SJ ' interesting applications in many areas of Physics. It is a (34- l)-dimcnsional Lorentz invariant field theory for a triplet 
t — of scalar fields n, living on the two-sphere 5'^ , = 1 , and defined by the Lagrangian density 

\o ■ 

^ ■ £ = Af 2 dnu ■ d^'n - 4 {dun A d^rif + ^ (d.n ■ d'^rif (1) 

' where the coupling constants and are dimcnsionless, and M has dimension of mass. The first two terms 
. correspond to the so-called Skyrme-Faddeev (SF) model [1], which was proposed long ago following Skyrme's idea 
• as the generalization to 3 + 1 dimensions of the CP^ model in 2 + 1 dimensions 3]. The interest in the SF 

. ^ I model has grown considerably in recent years since the first numerical knotted soliton solutions were constructed 
' ^, 5, 6, 7, 8J, as well as vortex type solutions [1, It has also found applications in Bose-Einstein condensates 
[ superconductors [H, [HI and in the Weinberg-Salam model In addition, it has been conjectured by Faddeev and 
Niemi [Ts'] that the SF model describes the low energy (strong coupling) regime of the pure SU (2) Yang-Mills theory. 
That was based on the so-called Cho- Faddeev- Niemi-Shabanov decomposition of the SU{2) Yang-Mills 

field Afj_, where its six physical degrees of freedom are encoded into a triplet of scalars n (n^ = 1), a massless U{1) 
gauge field, and two real scalar fields. The motivation for such decomposition originates from the fact that the triplet 
n can be used as an order parameter for a condensate of Wu-Yang monopoles, which classical solution can then be 
written as Aq = 0, Ai = difi A n, ? = 1, 2, 3, with n corresponding to the hedgehog configuration n — xjr. 

The conjecture of Faddeev and Niemi requires non-perturbative calculations to be proved or disproved and so the 
discussions in the literature have been quite controversial. Lattice field theory simulations have discouraged its validity 
[lO'l , and recently Faddeev himself has proposed some modifications for it [2^ . Gies [2l| has calculated the one- loop 
Wilsonian effective action for the SU (2) Yang-Mills theory, using the Cho-Faddeev-Niemi-Shabanov decomposition. 



'Electronic address: laf@ifsc.us p.br| 



t Electronic address: sawado@p h.noda.tus.ac.jp| 
■t Electronic address: kouichi@y ukawa.kyoto-u.ac.jp| 



2 



and found agreements with the conjecture provided the SF model is modified by additional quartic terms in derivatives 
of the n field. In fact, Gies obtained an effective action which up to first derivatives of n is of the form Similar 
results were obtained in [23] for the connection between the low energy Yang-Mills dynamics and modifications of 
Skyrme-Faddeev model using the Gaussian approximation of the vacuum wave functional. In this sense the extended 
version of the Skyrme-Faddeev model given by the theory ([T]) deserves some attention. Like the SF model it presents 
the internal 0(3) symmetry and admit solutions with topology given by the Hopf map ^ S"^. The main difference 
is the fact that ([T]) contains quartic terms in time derivatives and so its canonical Hamiltonian is not positive definite 
when hard modes of the n field are allowed. Indeed, that is a feature of many low energy effective theories. 

The theory ^ has many interesting aspects which deserves attention as we now explain. It has two sectors where 
the static energy density is positive definite and one of them intersects with an integrable sub-sector of ([T]) with an 
infinite number of conservation laws. In a Minkowiski space-time the static Hamiltonian associated to ([T]) is 



Hstatic = difi • 9^7? -I- ^ {difi A djfif 



(2) 



with i.j — 1, 2, 3. Therefore, it is positive definite for > 0, > and (3 < 0. That is the sector explored in [23j . 
and where static soliton solutions with non-trivial Hopf topological charges were constructed. 

One can now stereographic project S'^ on a plane and work with a complex scalar field u related to the triplet n by 



1 



(u + u* , —i {u — u*) , I M P — l) 



(3) 



It then follows that 



2i 



{dyud^u* - d^udyU*) 



Therefore, in terms of the field u the Lagrangian (jT]) reads 



(4) 



duU d^u* 



(1+ I n Pr 

The same static Hamiltonian ([2]) is now written as 

diU diU* 



{l+\u\^t 



(/3< 



il+\u\^f 



^static =4M^ 



{l+\u\^f ^ ' {l+\u\^f 



Therefore, it is positive definite for 



A/2 > : 



< 0; 



/3<0; 



j3e^ > 1 



(5) 



(6) 



(7) 



It is that sector that we shall consider in this paper. We show that it possesses static soliton solutions with non-trivial 
Hopf topological charge. If one compares ([T]) with the effective Lagrangian given by eq. (14) or (21) of [2l| one 
observes that, in order for them to agree, and /3 should indeed have the same sign. In the physical limit where 
the infrared cutoff goes to zero it follows that they both become positive, contrary to our assumption. However, the 
perturbative calculations in [2l| are not valid in that limit, and so not much can be said even about their relative 
signs. An interesting study of such sector was performed in [24] with an additional term in the action containing 
second derivatives of the n field and corresponding to the extra term obtained in ^21] when the hard modes of n are 
also integrated out. 

Another interesting aspect of the theory ([T|) has to do with its integrability properties. The Euler-Lagrange equations 
following from ([5]), or equivalently ((T|), reads 



{l+\u 

together with its complex conjugate, and where 



(8) 



JC^ = M d^u + — — 



(9) 
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The theory ([T]) has three conserved currents associated the 0(3) internal symmetry. However, using the techniques 
of [1^, [2^ one can show that the sub-sector defined by the constraint 

a^wa^u^O (10) 

possesses an infinite number of conserved currents given by 

where G is any function of u and w*, but not of their derivatives, and /C^^ is obtained from ^ by imposing pop . i.e. 

' e2 (1+ I u |2)2 

The fact that currents pT|) are conserved foUows from the identity Im (JC^d^^u*) — 0, the condition JC^d^u — 
foUowing from (fTU|) . and the equations of motion, which now read d^lC^ — 0. 

If in addition of pH]) one restricts to the sector where the couphng constants satisfy 

(3e^ = l (13) 

then the equations of motion simphfy to u = 0, and the theory becomes scahng invariant. That constitutes an 
integrable sub-sector of the theory ([1]) . Exact vortex solutions were constructed in [27| , using quite simple and direct 
methods. One can have multi-vortex solutions all lying in the same direction, and the vortices can be either static or 
have waves traveling along them with the speed of light. 

Apparently the integrable sub-sector defined by conditions (fTO)) and (fT3)) do not admit soliton solutions with non- 
trivial Hopf topological charges. In this paper we construct such solitons for the theory ([T]) in the range of the coupling 
constants defined in ([7]) and without imposing the constraint pop. We construct static solutions with axial symmetry 
by solving numerically the equations of motion. Due to the axial symmetry we have to solve two coupled non-linear 
partial differential equations in two dimensions. We choose to use the toroidal coordinates defined as 

1 ''0 r- 

X — y/z COS Lp 

P 

x^ — — \fz sin p — \ ~ COS ^ \J\ — z 



X 



P 

3 _ '^0 



— VT^sin^ (14) 



P 



where x*, i — 1,2,3, are the Cartesian coordinates in K'^, and (z, ^, (^) are the toroidal coordinates. We have 
^ z ^ 1, — TT ;^^^7r, 0^(^^ 27r, and ro is a free parameter with dimension of length. Notice that the usual 
toroidal coordinates have a coordinate ry > which is related to our z\>y z — tanh^ 77. 

We construct solutions which are invariant under the diagonal C/(l) subgroup of the direct product of the C/(l) 
group of rotations on the x^ x^ plane and the C/(l) group of internal phase transformations u — s- e*" u. We then use 
the ansatz 



y= /L5(£^g.e(.,?)+-^ (15) 

V 5(2,0 

with n being an integer, and < 51 < 1, — tt < < tt. 

By rescaling the Cartesian coordinates as — > x^ /tq, it is then clear that the equations of motion ([8]) will depend 
only upon two dimensionless parameters, namely (3 >\ and 

= -e^ rl > (16) 

since we are assuming ([T])- In addition, from ^ the static energy can be written as 

4M r. 



E = / d^xH^ 



+ 2(i?W + (/?e2-l)i?f) 



(17) 
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in terms of the dimensionless quantities 

E2 ^ M \ e\ I d^x 



(1+ I u 



.(1) _ _j_ f^3 j(g.-^r 1^ 

M\e \ J (1+ I u |2)4 



Ef' = / (18) 



^(2) _ 1 /"^a^ (a^.wa^.u*) 



^4 



M I e I 7 (1 



By replacing the ansatz ([15]) into the equations of motion ([8]) we get two non-hnear partial differential equations 
for the function 3(2,^) and 8(2;,^), depending upon the dimensionless parameters (3 and a. We solve those 
equations numerically usin g a standard relaxation method psj , with appropriate boundary conditions explained 
below. Derrick's argument [29| implies that the contribution to the energy from the quadractic and quartic terms 

must equal, i.e. E2 ~ 2 (^E^^ + (/3e^ — l) e'i^^ (see PT|) ) . That fixes the size of the solution and so the value 

of a. Therefore, in our numerical procedure we determine a on each step of the relaxation method using Derrick's 
argument. Our main results are the following. We find finite energy axially symmetric solutions with Hopf charges 
up to four. The axial symmetry comes from the use of the ansatz (jlSp . Therefore, our solutions are not necessarily 
the ones with the lowest possible energy for a given value of the Hopf charge. For the Skyrme-Faddeev model the 
solutions with Hopf charge 1 and 2 do present axial symmetry [1, d, [a 0i [1| ■ We believe the same may happen in the 
extended model ([T]). An interesting discovery we made is that a ^ as /3e^ ^ 1. That means that the solutions 
shrink in that limit, and disappear for /3e^ = 1. In addition, the energy also vanishes in the limit /Se^ —> 1. Since all 
three terms in the energy pT|) are positive they each vanish in that limit, including E^\ despite the factor (/3 — l) 

in front of it. Another interesting discovery is that E^^ is very small compared to the other two terms. In fact, for a 
wide range of the parameter we have 



Ei'^ 



E, 



+ 2 (i?« + (/3e2-l) £;f ) 
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Notice that Ej^'^ is a good measure of how close the solutions are of satisfying the constraint ((T0|) . Therefore our 
solutions are very close of belonging to the integrable sector possessing the infinite number of conserved currents 
given by dm)- It would be very interesting to perform calculations with a finer mesh to improve the precision of that 
ratio. One thing is clear however, the closer is fie^ to unity the better the constraint is satisfied by the solutions. 
However, the solutions shrink in that limit. The Hopfions then disappear and seem to give place to the vortex solutions 
constructed in [27] . That point certainly deserves further studies. It would be very desirable to perform simulations 
to see how the vortices of [27] behave when the conditions pUj) and are relaxed, and that fact is now under 
investigation. 

The paper is organized as follows: in section [H] we calculate the Hopf topological charges for the configurations 
within the ansatz l [T51) . in section UTTl we calculate a bound for the static energy (fT7)) using methods similar to those 
of references (32l . l33j] employed in case of the Skyrme-Faddeev model. The equations of motion in terms of the ansatz 
functions g {z,£^) and 8 {z,£^) are calculated in section [TVl and in section|V]we analyze the properties of the solutions 
in three limiting cases, at spatial infinity, at the circle + x| = r§, with X3 = 0, and also at the x^-axis. The 
numerical methods and solutions are described in section IVll 



II. THE HOPF TOPOLOGICAL CHARGE 

The static solutions of the theory ^ define maps from the three dimensional space to the target space 5"^. 
In order to have finite energy solutions we need the fields to go to a vacuum configuration at spatial infinity, and 
so we need n — > constant for r 00, with r being the distance to the origin of the coordinate system. Therefore, 
finite energy solutions map all points at spatial infinity to a fixed point of S^. Consequently, as far as the topological 
properties of such maps are concerned we can identify the points at infinity and consider the space to be instead 
of IR'^. Then the finite energy solutions define maps ^ S^, and those are classified into homotopy classes labeled 
by an integer Qjj called the Hopf index. Such index can be calculated through an integral formula as follows [13, |3l| . 



We first consider the mapping of the three dimensional space TR^ into a 3-sphere S^, parametrized by two complex 
numbers Zk, k = 1,2, such that ] Zi ]^ + ] 1^= 1. Using the ansatz (|15p we choose such map to be given by 

IR3^5|: Zi^yi^e^^ Z2 = V5e-'"'^ (19) 
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We then map such 3-sphere into the target S'^ by 



Zi 



(20) 



In fact, the mapping of the plane parametrized by the complex field u into the S"^, parametrized by the triplet of 
scalar fields n, is given by the the stereographic projection defined in ([3]). The Hopf index of such map is given by 

M 



with 



Qh 



1 

4^ 



k=l 



(21) 



(22) 



and where the integral is over the coordinates a;' of the three dimensional space TR^ . Even though A cannot be written 
locally in terms of u or n, its curvature can. Indeed, from (jH), (|20|) and ([22|) we have 



2 (l+|u|2) 



Therefore the Hopf index (PT|) can alternatively be written as 

1 



Qh - 



Performing the calculation using the toroidal coordinates (fT^ we find that 



A=P- 
ro 



-2 (1 - z) (1 - g) - ^^j=A d^O e^ + n^, 

VI — z V^ 



(23) 



(24) 



(25) 



where e^, and are the unit vectors in the direction of the changes of the position vector r under variations of 
coordinates, i.e. — dz?/ \ S^r |, and so on. The metric and volume element in toroidal coordinates (I14p are 



dz^ 



4z (1 - z) 



+ {l-z) dC^ + zdip^ 



d^x = 



1 fro 



2 \p 



dz d^ dip 



Therefore the Hopf index pT|) for the field configurations in the ansatz p5|) is given by 



Qh^^ C dz r d^ [dz (gd^Q) - {gdzQ)] 
2tt Jo 



We now impose the boundary conditions 

5(^-0,0 = g(z^ 1,0 = 1 for 

and 

Q {z,S. — — tt) = — m TT Q {z,^ = tt) = rmr for 

with m being an integer. It then follows that dzQ |^=-7r= dz& \^=-k= 0, and so from ((27| 

Qh = mn 



TT < < IT 



< z < 1 



(26) 

(27) 

(28) 
(29) 
(30) 
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III. THE ENERGY BOUND 



The static energy of the theory ([T]) in the regime ([7]) satisfy an energy bound very similar to that one found by 
Vakulenko and Kapitanskii :33] for the Skyrme-Faddeev model. The arguments are quite simple, and we start by the 
definition of Fij in (|23p which implies that 



2 _ ^ jd.udiu*) 



2 12 



and so 



diU diU* 



We now rewrite the static energy density ([6]) as 

diU diU* 



> 




AM 



,r2 '^1 



(l+|u|2 

Therefore, using (j32p and working in the regime ([7]) one has 



Hstatic > 4 W - _ (/3 - 1) i^: 



(31) 



(32) 



(33) 



(34) 



Notice that the bounds (|32p and ((34| are saturated by those configurations satisfying the constraint pO)l . Therefore 
we should expect the solutions to be driven in the direction of satisfying that constraint, and that is what we find in 
our numerical simulations. 

The static energy ([T7)) should then satisfy (assuming the regime (O) 



E > 2"/4j^^^iT3T (^l^^x^ 
We now use the Sobolev-type inequality [3^, [S^l 



1/2 



2M 
21/4 



d'xJF^^ 



2 - 1 / d-^xF, 



1/2' 



H 2 



d-'xjFfj 



(Fx Ft, > — 



1/1 



C V87r2 



where C is a universal constant. We then get the energy bound 

211/4 

where we have used (|24|) . The value of C found in [H, [s^l is 

1 = 83^/4^2^4 
o 

However, Ward in |33|] conjectures that a better value is 



.3/4 
iH 



3/2 



Taking Ward's value ([5^ one then gets the bound 



i;> 647r2 — V/3e2-l Q^/ 
e 



3/4 



(35) 



(36) 



(37) 



Notice that one should expect a decrease in the energy of the hopfions as (3e'^ approaches unity, and that is exactly 
what we find in our numerical calculations. 
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IV. THE EQUATIONS OF MOTION 



By replacing the ansatz (|15[) into the equations of motion ([8]) we get two partial non-linear differential equations 
for the functions g {z , ^) and 6 (z , ^). They are given by 



g{l - g) 4(z(l - z)f{pd,R, - R,d,p) + z{pd^R^ - R^d^p) 

r3 



(1 - 2g)p[-(4(z(l - z)YrA9 + zR^d^g) - g{l ~ g)(A{z{l - z)fSAQ + zS^d^Q + (1 - z)S^7i)'\ = (38) 



and 



7(1 - 5) [4(z(l - z)f{pd,S, ~ S^d^p) + z{pd^S^ - S^d^p) 
-3 



-(1 - 2g)p[-(4(z(l - z)fs^d^g + zS^d^g^ + .g(l - g)(4(z(l - z)Y RM + zR^D^Q + (1 - = (39) 

with p being defined in (|14p and where 



where 



and 



Rz 


:=-l[gil-g)z{l-z)^ 




+ l{va + Wfc)^ - ag{\ 


- g)vcdzQ 




-^(.9(1 -.9)^(1-^)4 


- a{Va - Vb) 


+ 7(^'a + Vbijd^g - ag{l 


- g)vcdiQ 


Rip 


:= —ag{l — g)vcn 






^avcd^g 


Sz 


= 5(1- 5)(5(1- 5)^(1 - 


z) - a{va ~ 


Vb) + l{Va + Vb)^dz<d + 


Si 


= g(l-g)(5(l- 5)^(1- 


z) - a{va - 


Wfc) + 7(wq + t^b))9^6 + 


^oiVcd^g 


Sip 


= "5(1 - 5) (5(1 - 5)^(1 


- z) - a{va 


- Vb) + 7(Wa + Wh)j 






a:=4^ 




7:=44(/3e2-l) 





(40) 



t'a :=i(4(z(l-z))2(a.g)2+z(ae9)') 

:= (.9(1 - g))2 (4(z(l - z))2(9,e)2 + ^(9^6)^ + (1 - z)n') 
vc := 5(1 - 5) (4(^(1 - z))\dzg){dze) + zid^g){d^e)) 



(41) 



Besides the boundary conditions (pS)) and leading to the Hopf index given by we impose as well the additional 
boundary conditions 



di9{z,0k=-^ = digiz,Ok=^ = 0, 
dze{z,0\z=o = dzQ{z,0\z=i = 0, 



for 
for 



< z < 1 
- vr < ^ < TT 



Notice that the equations of motion (|38p and (|39p are invariant under the transformations 

9iz,0^9iz,-0 e(z,o---e(z,-c) 

Therefore, we choose the boundary conditions 

%(2,0lc=o = o, e(z,$ = o) = o 



(42) 
(43) 



(44) 



(45) 



and perform our numerical calculations on the half-plane defined by < z < 1 and < ^ < tt. The functions on the 
other half-plane, namely < z < 1 and —n < £, < 0, are obtained from the symmetry 
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V. THE BEHAVIOR OF THE SOLUTIONS IN LIMITING CASES 



The solutions within the ansatz (jisp have axial symmetry around the a;'^-axis. In addition, from ([3]) and (jlSp one 
observes that the condition 713 = constant implies g = constant. The numerical solutions satisfying the boundary 
conditions ([^51) . and (I45p . that we find are such that for a given value of ^, 5 is a monotonically 

increasing function of z. In addition, for a given value of z, g does not vary much under variations of ^. It then turns 
out that the surfaces in IR'^ corresponding to constant have a toroidal shape. They are in fact deformations of the 
toroidal surfaces obtained by fixing z and varying ^ and ip (see (|14|) ). It is then important to analyze the behavior of 
the solutions in three regions, namely the spatial infinity and close to the a;'^-axis where due to the boundary condition 
(|28p we have — 1, and close to the circle on the a;^-plane of radius ro and centered at the origin, where = —1. 
We now perform such analysis. 



A. The behavior at spatial infinity 



From (fT4l) one observes that 



7^ ^ Xi ~r ~r 



(1 + z COS^) 



(i-VT 



^0 



Therefore, in the toroidal coordinates (jT4|) the limit r ^ 00 corresponds to the double limit z and f 0. We 
choose to perform such double limit through the parametrization 



z ~ £1 cos <7 



^ ^ Si sin cr 



with 



£1^0 0<cr<- 

- - 2 



(46) 



That corresponds to polar coordinates on the plane (y^, ^), and the derivatives are related by 

1 



cos cr 



(e cos crde — sin a da) 



Sm (T H Oct 



Taking into account the boundary conditions ([28| . (|43ll and ([45]) we then assume the fields to behave as 
where si and ri are constants, and with the test functions satisfying 



/(f) 







/'(o)-o 



h{Q) = Q 



'■(f) 



= 



(47) 



(48) 



where primes denote derivatives w.r.t. a. Replacing the expansions (|47)) into the equations of motion ((38|) and (|39p . 
and keeping only the leading terms, i.e. those with the smallest possible power of £1, we get the following equations 
for the test functions 



1 (f 



r , 1 



and 



) + o ( y ) - tancr — + - si (s 
/' 



2n' 



= 



/ 



tancr 



h' 



+ ri in + si - 1) = 



A solution of satisfying the boundary conditions (pS)) is given by 

/ - (cos cr)^'"' with Si = 2 (I n I +1) 

Replacing ([5T|) into ([50)) and performing the change of variables 



(49) 



(50) 



(51) 



X = sm a 



< a < 



< a; < 1 
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we get that the equation (|50p becomes an hypergeometric equation of the form 

x{l~x) dlh + i [1 - (2 I 71 I +3) x] dji + ^ ri (ri + 2 | n | +1) = (52) 



The boundary conditions (|48p imply that 

h \x=Q= dxh \x=i^ finite (53) 

There are two independent solutions for equation (j52|) . namely 

h2 = x'/'F(\n\+'-l + l,-'-l + l,l,x 

where F {a, f3,j,x) is the Gauss hypergeometric series. One can check that such series does not converge since in 
both cases we have a + /3 — 7=| n |> 1. We then have to choose ri in order to truncate the series. That leads us to 
the allowed values for ri 

ri^2l + l / = 0,1,2,... 

Since the first term of F is 1, we see that only /12 satisfies the boundary condition (|55|) . Therefore, the allowed 
solutions for O are 

e - £f+^ sinai^ (^1 n I +/ + ^,-/,^,sin2CT^ Z = 0,1,2,... 

In our numerical calculations only solutions corresponding to / = appear (so F — 1), and therefore the behavior of 
our solutions at spatial infinity, i.e. z , ^ ^ 0, is given by 

5~£j'l"l+') (cosa)'l"l~zl"l(z + e') e-eisind-e (54) 

That implies that the u field at spatial infinity behaves as 

1 y,2|n| + l 



(2ro)l"l+' 

with p'^ = Xi + X2, and r"^ = x\ + x\ + x\. Therefore, we indeed have n (0, 0, 1), as r ^ 00. 



(55) 



B. The behavior around the circle x\-\- X2 ~ r^, and x-i — Q 

The circle x\ + x"^ — ^O' ~ ^ correspond to z = 1 and any ^ and ip, and so we write 

£2 = 1 - z -> 

We assume the following behavior for the fields 

9^l-e^^F{0 Q^HiO-e^KiO (56) 

with S2 and r2 being constants. The boundary conditions ([29|) . (|42|) . ((43|) and (|45ll impose the following conditions 
on the trial functions 

d^g |5=0,±7r= ^ F' |^=0,±7r= 

e 1^=0= ^ H |^^o= K \^=o= 

a.e u=i= ^ r2 > 1 (57) 
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The last condition plays an important role in our analysis. In addition, we do not want the z-derivative of g to diverge 
close to z = 1, and so we need S2 > 1. Replacing those expansions in the equations of motion ([55)1 and ([55]) and 
keeping only the leading terms we get 



1 /F'\' 1 f F' 



2 



F 



2 ~F 



F 



- 4 {H'f 

1 fF'\' 1 fF 



(58) 



and 



H" F' ^ 



aF' 



l + iH'Y 



F 



1 fF' 



F 



1 - {H'f G_ + a [F' {H' 



F 



~F^2~Flr 



1 F' H" 
G'_ + --(G+ + G-) + — G- 



(59) 



where the primes denote derivatives w.r.t. ^, and where (5^2,1 is the Kronecker delta and so those terms only exist in 
the case S2 = 1- In addition we have introduced the quantities 



G± EE F 



1 fF' 



F 



1 fF' 



and a and 7 are the values of the functions (|40p in the limit z — s- 1, i.e. 



F 



1 



Notice that the function K (^) introduced in ()56|1 is not determined in such expansion up to leading order. 

For the numerical calculations we perform, the relevant solution for equations (|58p and (|59p . satisfying the boundary 
conditions (l57l) , is 



S2 =\ m 



F — const. 



and so 



g ~ 1 — const. (1 — z 



\m\ 



H = mi 



^ 



Therefore, the behavior of the u field around the circle x\+ x\ = r^, and = is 

M const. (l-z)l"l/^ e'('"«+"'^) z^l 
Therefore we indeed have n (0, 0,-1) on that circle. 



(60) 



C. The behavior of the solutions around the x'^-axis 

^From (|14p we see that the x'^-axis corresponds to z = 0, and any ^ and (p. Therefore we denote 

£3 = z ^ 

We assume the following behavior for the fields 

g^efJ^iO e^mO + el'lCiO (61) 



with S3 and being constants. The boundary conditions ([29]), (|42|). ([43]) and (|45)) impose the following conditions 
on the trial functions 

lc=0.±7r= ^ T' |^^o,±7r= 

e 1^=0= ^ n |e=o= o /c 1^=0= 

5.9 U=o= ^ r3 > 1 



11 



The last condition will play an important role in our analysis. In addition, we do not want the ^-derivative of g to 
diverge close to z = 0, and so we need S3 > 1. Replacing those expansions into the equation of motion psp and 
keeping only the leading terms we get 

(n^ - si) [1 + Ss.sJ'l {4 + - [4 + n") 5s,,iTa {si - n^) = (62) 

where a and 7 are the limiting values of the functions (|40p in this expansion, i.e. 

a = 47(l-cosO' 7= 4 (/5e2-l) (l-cos^)' 

So, apparently the only reasonable solution for ([62]) is 

S3 =1 I 

Now replacing the same expansions in the equation of motion ([55)) for the cases S3 > ra > 1 and 1 < S3 < rs we get 
the equation 

-T_/// '77/ O 

^ + ^- — = S3>r3>l l<S3<r3 (63) 

n, J' p 

and so the solution is given by •^J^ = const., where p is the liming value of the function introduced in (jl4p . i.e. 

p —1 ~ cos ^ 

For the case 1 = S3 < rs the equation of motion ([55]) gives instead 

(I + 27.F) — + (l + 4(a + 7).F) — -(I-27.F) l = S3<r3 (64) 

rt J- p 

We then notice that even though the function K, introduced in ([61]) does not enter into the equations, the exponent 
r-i related to it does play a role in the expansion up to leading order. The difficulty of this case however is that the 
expansion up to leading order is not enough to determine the trial functions T and 7i separately. We would have to 
go to the next to leading order terms to get those functions. 



VI. THE NUMERICAL ANALYSIS 



The equations ([351) and with the boundary conditions (US]), (HI]), (gH), (gS]) and (gS]) can be solved numerically 
by the standard successive over-relaxation (SOR) method [28*1. Such well known method is a powerful tool for finding 
solutions of this kind of complex set of equations. In order to find a solution / of an elliptic equation = p where 
£ represents some elliptic operator and p is the source, we rewrite it as a diffusion equation, 

^=uj{Lu-p). (65) 

The idea is that an initial distribution u relaxes to an equilibrium solution / as t — s- 00, and w is called the over- 
relaxation parameter. (Normally one chooses 1 < < 2 for faster convergence). We apply that algorithm to our 
equations by putting p8|) . p9|) into right-hand side of the diffusion equations ((65)) . One thing we should care about 
is the sign of the elliptic operator. In order to get convergent solutions we found that we must choose the r.h.s. of 
()65p in such a way that the second derivatives of g and 0, w.r.t. z and ^, appear with a positive sign. Therefore, we 
have used the diffusion equations as 

- = -0;^ —^u^B (66) 

where A and B stand for the expressions on the l.h.s. of (|55)) and (|55)) respectively, and w was taken to be unity in 
most of the simulations. 

In order to estimate the derivatives in the equations, we have discretized the equations using central differences on 
a rectangular cubic lattice (z,^). We have investigated the cases of various mesh sizes and found a good convergence 
for a mesh size (iVz, = (80, 80), i.e. we divided the intervals < z < 1 and < ^ < tt into 80 segments each. 
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Qh 


(m,n) 


a 


E 


E/ -Ebound 


1 


(1,1) 


0.53 


236.2 


1.18 


2 


(1,2) 


0.79 


390.6 


1.16 




(2,1) 


0.66 


466.7 


1.39 


3 


(1,3) 


1.09 


548.7 


1.21 




(3,1) 


0.80 


715.7 


1.57 


4 


(1,4) 


1.41 


697.3 


1.23 




(2,2) 


0.94 


765.7 


1.36 




(4,1) 


0.99 


965.7 


1.71 



TABLE L The static energy E given by (|17p (in units of -r—r) s-nd lli6 parameter a, defined in (|16|l . for /3e^ — 1.1. We also 
compute the ratio of the energy E to the bound given in (I37[) . i.e. Ebound — 64%^ ^ \/f3e^ ~ 1 Q^^. 



The parameter a introduced in (|16p is determined by using Derrick's scaling argument. Indeed, that argument 
implies that the contribution to the energy coming from the quadratic and quartic terms in derivatives should equal, 



I.e. 



from pT]) one has E2 = 2 (jS^^^ + i^fie^ — l) E^^^^ . Therefore, we choose an initial value for a, and on each step 



of the relaxation method we determine a new value for a by imposing that condition. By using the ansatz 
finds that the energy pT|) can be written as 



E 



AM 



dz 



L a\ 



3 one 



(67) 



with a being defined in (jl6[) and 



£2(^,0 = TT 



.(1) 



W(l - 5)z(l - z) 

P{{va - Vbf + vl) 



(5(1-5))^(^(1-^))^ 
{g{l-9)Y{z{l-z)Y 



(68) 



where Wa , Vb and are defined in (|4ip . Consequently, the value of a on each step of the simulation is calculated by 
the formula 



a =2 



dz d^ (e« (z, + iPe' - l)ef ^ (z, 0) 



IodzJ^^d^e2iz,^) 



(69) 



We have made several simulations studying how the solutions depend upon the integers n and m (introduced in the 
ansatz (|15p and boundary conditions (j29p respectively) and the parameter /3e^. For the particular case of /3e^ = 1.1, 
we have calculated numerical solutions with Hopf topological charge up to 4. The energies as well as the size of the 
solutions, measured by the parameter a, are given on the TableHl We also compare the energies with the bound given 
by ((37|) . Very probably the solutions with Hopf charges 3 and 4 correspond to excited states. For the Skyrme-Faddeev 
model the minimal energy solutions with Hopf charge Qh > 2 exhibit quite complicated structures like, twists, linking 
loops and knots, than the simple planar tori [5j, Q. Thus, the solutions with the minimum value of energy on those 
sectors very probably do not have axial symmetry. The solutions corresponding to (m, n) = (1, 1) and (m, n) = (1, 2), 
and so Hopf charges 1 and 2 respectively, may correspond to the minimum of energy. That is in fact what happens 
in the Skyrme-Faddeev model where the solutions of charge 1 and 2 do present axial symmetry. 

We have made a more detailed analysis of the solutions corresponding to (m, n) = (1, 1) and (m, n) = (1, 2). We 
have studied how those solutions depend upon the parameter /?e^, and the results are quite interesting. The energies 
of those two solutions decreases as /?e^ approaches unity from above, and they vanish at /3e^ = 1. The same is true 
for the size of the solution. The parameter a decreases as /?e^ ^ 1 and vanishes at /3e^ — 1. That means that the 
solutions cease to exist at that point. The results are shown in the plots of Figure [1] In addition, the contribution to 
the total energy from the term E^'' (see (IT71) and (HH)) is very small for the range considered for the parameter /3e^. 
It corresponds in fact to about 0.1% ~ 0.5% of the total energy (see Figured]). From ^TE^ we see that the term iJ^^'' 
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is a good measure of how close the sohition is of satisfying the constraint (|10p . which leads to the infinite number of 
conserved currents (|lip . We see that those two solutions are very close of belonging to that integrable sector. The 
shrinking of the solutions can perhaps be understood by Derrick's argument. For some reason which is not very clear 
yet the dynamics of the theory keeps the quartic term of the energy iJ^^^ very small. We have pointed out in section 
mil that configurations that satisfy the constraint (|10p have energies closer to the bound p7|) . and that may be a way 
of understanding why £'4^' is so small. The other quartic term, namely £^4^' is multiplied by (/3e^ — l) (see ([T7)) ) 
and so its contribution vanishes as /?e^ — > 1. Therefore, we are left with the quadratic term E2 only, which cannot 
be balanced by the quartic terms anymore. Since that term scales as £2 — > A £2, as a;* — *■ Xx^, the solution tends 
to shrink. Another fact is that the equations of motion in the sector of the theory where /?e^ — 1, and where (fTU]) 
is satisfied, present scale invariance (see comments below ([7])). Therefore, if a localized solutions like the hopfions 
exists, then any re-scaling of it would also be a solution. But that is very improbable to happen in such a theory. 
Consequently, only vortex solutions like the ones constructed in [S^l may exist in that sub-sector. 

The solutions for the ansatz functions g {z,£^) and 8 {z,£,) satisfying the boundary conditions (|28p . ([29| . (|42ll . (|43|) 
and (|45ll are quite regular. In Figures [2] and [3] we give the plots of those functions for the cases {m,n) = (1, 1) and 
{m,n) = (1,2) respectively, and for the values /3e^ = 1.1 and 5.0 in both cases. We show the plot in the range 
< C ^ ■"'j s-iid the solution in the range — tt < ^ < 0, is obtained through the symmetry ([44]) . Notice that the 
changes of g {z, ^) and Q (z,^) under variations of /? are very small. In order to visualize those changes better we 
also present in Figure [3] cuts of the function g at the borders ^ = and £, = tt, and for several values of (3 e^. We also 
show in Figure [5] cuts of Q at the borders z = and z = 1 for several values of /3 . 

We also present the densities of the static energy ([T7]) in the cylindrical coordinate space (p := a^/z/p, x"^) in Figl6] 
for Qff = 1, and in Fig. [7] for Qh = 2, for several values of /3e^. The following relations between two coordinates are 
quite useful to visualize them: 



a2 _ p2 _ (^3)2 



2^2 



L4 - 2a2(p2 _ (2.3)2) ^_ (^2 + (^3)2)2 J ' ~ (^2 + ^2 ^ (^3)2)2 " ("0) 

The energy density of (m, n) = (1, 1) {Qh ~ 1) solution has the lump shaped, and of the [m,n) ~ (1,2) {Qh = 2) 
solution exhibits the toroidal configuration. For both case, one easily observes the shrinking of the solutions as 
/3e^ 1, and so an increase of the density around the origin. 

In Fig[51 we display the energy densities for solutions with charges Qh > 2. The n = 3,4 solitons exhibit the 
toroidal shape. The radius of the tori increases as n grows. On the other hand, all n = 1 solitons seem to be lump 
shaped, however, for m ^ 2 they have a depression close to the origin. The size grows especially in the x^ direction 
as the charge increases, and as that happens a second peak gradually emerges. 
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gH = l, [(m,n) = (l,l)] Qh = 2, [(m,n) = (1,2)] 

FIG. 1: The plots on the Lh.s. correspond to the solution where (m,n) = (1, 1) and so Hopf charge equal to 1, and those on the 
r.h.s to the solution where (m,n) = (1, 2) and so Hopf charge equal to 2. On the first and second rows we plot the total energy 
E (see (|17p ) and the parameter a (see (|16|) ') respectively, against (3e^. On the third and fourth rows we show the fraction of 
the total energy corresponding to the two quartic terms, E\ and El respectively, as a function of /3 e . Notice that due to 
Derrick's scaling argument we must have E2 = 2 ^£'4^'' + (/3 — l) £4^'^ , and so the sum of the values shown on the plots on 
the third and fourth rows gives 0.5. 
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FIG. 2: The ansatz functions g{z,(,) and e(«,{) for the solution (m,n) = (1, 1) {Qh = 1) for /3e^ = 1.1 (left), and for /3e^ = 5.0 
(right). 
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FIG. 3: The ansatz functions g{z, and e{z, for the solution (m, n) = (1, 2) {Qh = 2) for /3e^ = 1.1 (left), and for /3e^ = 5.0 
(right). 
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(a) g{z,0 at ^ = (left) and ^ = tt (right); 




(b) 6(^,0 at z = (left) and z = 1 (right); 

FIG. 4: Cuts of the ansatz functions g(z,^) and Q{z,^) at the boundaries for the solution (m, n) = (1, 1) {Qh = 1) and for 
some values of /3 shown on the plots. 
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(a) g{z,0 at ^ = (left) and ^ = tt (right); 




(b) 6(^,0 at z = (left) and z = 1 (right); 

FIG. 5: Cuts of the ansatz fiinctions g{z,^) and Q{z,^) at the boundaries for the solution (m, n) = (1,2) {Qh = 2) and for 
some values of /3 shown on the plots. 
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FIG. 6: The energy density (in units of — ) for the solution (m, n) = (1,1) {Qh = 1) on the cylindrical coordinates plane 
(p := a^/z/p, x^). 
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FIG. 7: The energy density (in units of — ) for the solution {m,n) — (1,2) (Qh ~ 2) on the cylindrical coordinates plane 

1^1 

(p — ay^/p,x^). 
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(2,2) (2,1) 




